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Abstract 

Magnetic insulators can be characterized by a gap separating the singlet ground state from 
the lowest energy triplet, S = 1 excitation. If the gap can be closed by the Zeeman interac- 
tion in applied magnetic field, the resulting S = 1 quasiparticles, triplons, can have concen- 
trations sufficient to undergo the Bose-Einstein condensates transition. We consider macro- 
scopic properties of the triplon Bose-Einstein condensates in the Hartree-Fock-Bogoliubov 
approximation taking into account the anomalous averages. We prove that these averages 
play the qualitative role in the condensate properties. As a result, we show that with the 
increase in the external magnetic field at a given temperature, the condensate demonstrates 
an instability related to the appearance of nonzero phonon damping and a change in the 
characteristic dependence of the speed of sound on the magnetic field. The calculated mag- 
netic susceptibility diverges when the external magnetic field approaches this instability 
threshold, providing a tool for the experimental verification of this approach. 
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1. Introduction 

Macroscopic systems governed by quantum mechanics of interacting particles, with en- 
sembles of cold atoms and quantum magnets being the intensively investigated examples, 
attract a great deal of interest. As well understood, they have interesting similarities, show- 
ing the common physics of these two seemingly different realizations. Macroscopic ensembles 
of bosonic atoms can undergo Bose-Einstein condensation (BEC). The properties of these 
atomic condensates in different systems became one of the most interesting fields in modern 
experimental and theoretical physics. On the other hand, triplet quasiparticles in magnetic 
insulators, being bosons with spin equal to one, can undergo the BEC transition. The first 
example of this kind of condensate is the systems far away from the equilibrium, where 
high concentration of spin excitations (magnons) is achieved by a strong resonant external 
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Figure 1: (Color online) Evolution of the magnon states with the increase in the magnetic field. When the 
bottom of 5*2 = — f band has negative energy, magnons (left panel) become "triplons" (right panel). 



magnetic field pumping [iHJ]. At a sufficiently intense pumping, this type of condensation 
can occur at temperatures of the order of 100 K. Second example j5| is the magnons in 
superfluid He^, where the magnetism appears due to the very small nuclear rather than 
electron magnetic moment. This process occurs at temperatures lower than 10~^ K. Third 
interesting example is presented by triplons, appearing if the gap in the triplet excitations 
spectrum separating it from the singlet ground state, can be closed by the Zeeman effect 
of the magnetic field [gI, 0] • The field splits the spin S = 1 excitations into three branches 
with Sz = 0,±1. When the Zeeman shift of the 5*2 = —1 branch exceeds the initial gap 
in the spectrum, a finite population of triplons, which can be considered as the reconstruc- 
tion of the ground state, is formed even at zero temperature, as shown in Fig{TJ The BEC 
transition is seen as the transverse magnetization, which occurs if the field becomes higher 
than the critical one jsj. The Bose-Einstein condensation of triplons was experimentally 
first observed jsl, |lo[ and thoroughly studied in the magnetic insulator TlCuCla jioj where 
magnetic properties are due to the presence of noncompensated single-electron spins of Cu 



ions. Similar findings in a variety of other compounds followed shortly |lll-ll5|. Recently, 



Bose-Einstein condensation of magnetic excitations was studied theoretically in compounds, 
where the magnetism is due to the spins of vanadium 16| or chromium 17|] ions. These 
uniform, three dimensional with various degrees of anisotropy and easily tunable systems 
provide the researchers with new abilities to study the BEC, including the effects of disor- 
der. The tunability is realized by easily changing the chemical potential with the external 
magnetic field H^xt- The condensates were studied in a variety of magnetic fields up to the 
fields where the Zeeman effect strongly changes the properties of the systems due to the 
magnetization saturation. 

Since the triplon BEC occurs in solids, it possesses at least three interesting features 
resulting from its coupling to the host lattice. First, it is directly influenced by the spin- 



orbit coupling inherent to the solid where it is located [18|. For atomic BEC the analog 



of spin-orbit coupling can be produced by special combination of laser fields in some cases 
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Second, the triplon BEC is coupled to phonons, and, therefore, can provide a test for 



the effects of decoherence and hysteresis due to the couphng to the lattice [20|, |21|. Third, 
an exter nal p ressure applied to the crystal can strongly change the properties of the triplon 
BEC 2^,l23|. In addition, the solid-state background for the triplon BEC makes it accessible 
by a variety of experimental techniques, not applicable for the atomic condensates. 

Theoretical description of triplon BEC in relatively weak fields, where the Zeeman split- 
ting is much less than the width of the triplon band in the Brillouin zone, and, correspond- 
ingly, concentration of triplons is small, can be done in terms of weakly interacting Bose 
gas theory. The observation of the sound-like Bogoliubov spin-flip mode, the fingerprint of 
the BEC in an interacting Bose gas, in the neutron scattering experiment 2^ provided the 
strong support of this approach to the BEC of triplons. 

Main tools for these studies are based on the approach usually referred to as the Hartree- 
Fock-Popov (HFP) approximation, neglecting the anomalous density term^. The drawback 
of this approximation is the prediction of a jump in the triplon density, and, in turn, in the 
sample magnetization across the transition temperature. Although this approach provides 
a good quantitative description of the specific heat C„, it is unclear whether it captures 
qualitatively all relevant physics of the condensate. For this reason it will be interesting 
to study the behavior of the triplon BEC using a more accurate approach in order to see 
whether the analysis beyond the HFP approximation reveals new qualitative features. 



In our previous work (25[ we have shown that an extended version of the mean field 



approach (MFA) taking into account anomalous density terms in the condensate, improves 



the situation considerably [26|. The key part of this approach concentrates on finding the 
speed of sound-like Bogoliubov excitations in the BEC of interacting particles. This Hartree- 
Fock-Bogoliubov (HFB) approximation 27j leads to a continuous magnetization across the 
transition, in agreement with the experiment. By applying this approximation to the BEC 
of triplons we have shown that when the external magnetic field i^ext exceeds a critical value, 
H^^^, the speed of sound becomes complex and the BEC state undergoes an instability. 

Another general feature clearly seen in the triplon BEC is the dependence of its physics 
on the bare magnon dispersion ^k- The non-parabolic bare dispersion of magnons 2l|, |28 - 



30| leads to a non-trivial dependence of the transition temperature T^. on the concentration 
p ~ M(T, i^ext) and, hence, on -ffext- The bare dispersion, being itself i^ext-independent, 
determines the interplay of kinetic and potential energy of a macroscopic system, and, 
therefore, plays a crucial role in the BEC properties. The effects of the bare dispersion are 
clearly seen experimentally in TlCuCla as the p-dependence Tc ~ p'^^'^K The exponent 0(p) 
approaches 2/3 at low concentrations (low Tc), jilj as expected for the parabolic Ek ~ k"^, 
while at T > 2.5 K, (f){p) is close to 0.5. We will address the role of the dispersion below in 
the paper and present simplified for the parabolic dispersion calculations. 

In this paper we address the macroscopic properties of the BEC in terms of the speed of 
sound, specific heat, and magnetic susceptibility, in the HFB approximation and show both 



^Although connecting this approximation with the name of Popov in the hterature is not adequate since 
Popov never neglected cr-terms (V.N. Popov, Functional Integrals in Quantum Field Theory and Statistical 
Physics, Reidel, Dordrecht, 1983) we will use the acronym HFP approximation for the historical reasons. 
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the quantitative and qualitative importance of the anomalous averages. This paper is orga- 
nized in the following way. In Section II we outline main features of the HFB approximation. 
In Section III we discuss the specific heat and magnetic susceptibility in this approximation 
and show the importance of the anomalous averages. In Section IV we provide a detailed 
analysis of the instability caused by external magnetic field in terms of the speed of sound 
and magnetization. This instability is the qualitative effect, not expected if the anomalous 
averages are not taken into account. The Conclusions summarize our results. 

2. The Hartree-Fock-Bogoliubov approximation. 

We begin with the Hamiltonian of triplons as a nonideal Bose gas with contact repulsive 
interaction: 

V^t(r)(K-;x)V'(r) + |(vt(r)V^(r))' , (1) 

where ip{r) is the bosonic field operator, U is the interaction strength, and K is the kinetic 
energy operator which defines the bare triplon dispersion in momentum space. Since 
the triplon EEC occurs in solids, we perform integration over the unit cell of the crystal 
with the corresponding momenta defined in the first Brillouin zone. Below the spectrum 
will be assumed isotropic: = e{k). The parameter fi characterizes an additional direct 
contribution to the triplon energy due to the external field and being rewritten as 

/i = gfisHc^t - Ast (2) 

can be interpreted as a chemical potential of the Sz = —I triplons jsj. In Eq. g is 
the electron Lande factor, /xb is the Bohr magneton and Ast is the spin gap separating the 
singlet ground state from the lowest-energy triplet excitation (Fig{T]). 

In general the Hamiltonian in ([T]) is invariant under global U{1) gauge transformation 

'0(r) ^ e^"V'(r) (3) 

with a being a real number. However, this symmetry is broken in the condensed phase, 
where T < Tc, and restored for the normal phase, T > Tc. Note that in the experimental 
studies of BEC of atomic gases the density p (or equivalently, the total number of atoms) 
is fixed initially and the chemical potential n should be calculated self-consistently, while in 
the case of triplon BEC the chemical potential is fixed by the magnetic field and the density 
should be determined within an appropriate approximation. 

Here and below we adopt the units /c^ = 1 for the Boltzmann constant, ^ = 1 for 
the Planck constant, and V = 1 for the unit cell volume. In these units the energies are 
measured in Kelvin, the mass m is expressed in K~^, the magnetic susceptibility x for the 
magnetic fields measured in Tesla has the units of K/T^, while the momentum and specific 
heat Cy are dimensionless. Particularly, the Bohr magneton is hb = he/2mQC = 0.671668 
K/T, where mg is the free electron mass, and e is the fundamental charge. 
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2.1. Condensate phase T <Tc 

As it has been shown by Yukalov 27| the spontaneous gauge symmetry breaking is the 
necessary and sufficient condition for Bose-Einstein condensation and can be reahzed by the 
Bogohubov shift in the field operator as 

ij{r)=v{r)+i^{r), (4) 

where for the uniform case the function v{r) is a constant defining the density of condensed 
particles as 

Po = (5) 
Since by the definition of the average of ip'^{r)ip{r) is the total number of particles: 

N = I (iV(V^^(r)^(r)). (6) 



with the density of triplons per unit cell p = N/V, from normalization condition 

P = Po + Pi (7) 

one immediately obtains 



^ j^d\{^\v)ij{v)). (8) 



V 

Therefore the field operator ip{v) determines the density of uncondensed particles. Note 
that ?/^(r) and the condensate function should be considered as independent variables being 
orthogonal to each other in terms of the expectation values: 

dVV'(r)?;(r) = 0. (9) 

Consequently, passing to the momentum space 

i,{v) = ake^'^r ^ j "ke^''^ (10) 

and inserting (fTOl) into ([8]) one can easily see that in (fTOi) the summation by momentum in 
the finite volume should not include k = states. Below we indicate this rule by introducing 
prime sign in the momentum summation. 

Now using (jlj) and ffTOl) in ([1]) we present the Hamiltonian as the sum of five terms 



H = Y.Hn, ill] 



n=0 

labeled according to their order with respect to and a|^. The zero-order term does not 
contain the field operators of uncondensed triplons 

^0 = ^Po - (12) 
^ 5 



The first order term is identically zero 



^1=0, 

due to the orthogonality condition ([8]). The second-order term is 



H2 

The third-order 



{6k - n + 2U po) alai. 



(13) 



(14) 



(15) 



as well as fourth-order 



k,p 



u 



(16) 



k,p,q 



terms are rather complicated and a diagonalization procedure is needed. For this purpose 



in HFB approximation the following procedure is usually implemented [32 



(17) 
(18) 



where (aj^Op) = ^k.p'^k, (a-k^p) = 5k,-pC"k with rik and being related to the normal (pi) 
and anomalous (a) densities as 



Pi 



k 



k 



(19) 



Here we underline that the main difference between the HFP and HFB approximations 
concerns the anomalous density: neglecting a as well as (flkCtp) in fll4p . fll7p . and (1181) one 
arrives at the HFP approximation, which can also be obtained in variational perturbation 
theory 33|]. However, the normal, pi, and anomalous averages, a, are equally important 
and neither of them can be neglected without making the theory not self-consistent 341-136 . 
Altough pi and a are functions of temperature and external magnetic field, we omit explicit 
dependences in the formulas when it does not cause a confusion. 

Taking into account that = (H^) = 0, the expectation value of the energy can be 
evaluated as E(T < Tc) = Eq + E2, where 

Up'o 



En 



U 



PPo 



and Eo 



2 " 2 

{H2), with H2 being quadratic in a^, aj^: 

A 



Ho 



E 



— [ flkfl-k 
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a 



4a-k 



(20) 



(21) 



In the last equation 



(22) 



A = U{po + a). (23) 
The next step, in both approximations, is the Bogohubov transformation 

Ok = Mk^k + ^'k&^k' «k = "k^k + ^'k&-k (24) 

to diagonahze H2. The operators 6k and can be interpreted as annihilation and creation 
operators of phonons with following properties: 

K bl] = 5k,p, (&L&-k) = (&k&-k) = 0, (25) 
(bib^) = fsiEu) = (26) 

where (3 = 1/T. To determine the phonon dispersion Ek we insert ( |24|) into (12T|) and require 
that the coefficient of the term b^b-k + ^l^k^k vanishes, i.e: 

(^kUkVk + Y (wk + ^k) = 0- (27) 
Now using the condition ul, — v^ = l and presenting u^, as 

yields 



Jul - El = -A (29) 
that is 

Ek = ^{uJk + A){uJk- A) (30) 
where and A are given in Eqs. (122]) and 



Further requirement for Ek follows from the Hugenholtz-Pines theorem [37[: at small 
momentum k the spectrum should be gapless, and, therefore, the phonon dispersion is 
linear: Ek ^ ck + 0{k'^), where c can be considered as the sound speecil. This linearity can 
be achieved by setting 

cjfc - A = £k, (31) 



^It can be shown that 3J, |38[ A is related to the normal (En) and anomalous (Sa) self-energies as 
En = A + /i and Ea = A, respect ive l y. T he last two equations give En — Ea = jjL which is again in agreement 
with Hugcnholtz-Pines theorem 
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which together with Eq. (l23|) yields 

^les = ^^- 2Up = -A. (32) 

With this choice one obtains 

Ek = V^V^k + 2A, (33) 

with the sound speed 

c=./5, (34) 
V m 

i.e. A = mc^. Here m has the meaning of the triplon effective mass, characterizing the 
dispersion in the hmit of small momenta Ek ~ k"^ /2m. 

2.2. Condensed fraction and the condensate energy. 

Having fixed and one can find normal and anomalous densities as well as the energy 
by using equations ffTU]) -f l2S]) . For example, inserting into f lT^ results in 

k k k 

where we used the relation Ukt'k = —^/2Ek and introduced notation Wk = 1/2 + fsiEk). 
Similarly one obtains: 

Pi = Ei-W) = E - 1) • P'') 

S = (ft) = Y. EMEt) + lY(Ek- uit), (37) 

k k 

with Uk = Ek + A. 

When the bare dispersion Ek is isotropic, the momentum summation can be done with: 

k (271") ^0 

As a result, at T = the quantities in Eqs. fl5^ -f lTr|) can be represented as: 

(40) 

iEk + A)k^dk], (41) 



Pi(0) = 


-E 

k 


V Ek 


a(0) = 


A 


1"^ k'^dk 


47r2 ^ 


lo Ek 


^2(0) = 


■ 47r2 ( 


/ Ekk^dk 
Jo 
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The divergences in these integrals can be regularized by introducing a cutting parameter A 
(A — oo at the end of the calculations) or equivalently by using the dimensional regular- 
ization scheme. Now, assuming for the moment that, for T = case Sk = k"^ /2m and using 
dimensional regularization one obtains: [f| 



Pi(0) = (42) 



(Am)=^/2 
3^ 

a(0) = 3pi(0), (43) 
Summarizing this subsection we rewrite the above formulas as: 

Pi = o,2 + / 7:r^fB{Ek)—^ — , (45) 



3vr2 J (27iy'^ ' E, 



k 



a ^ - A J ^ME,)-, (46) 



where Bose distribution of phonons {Ek) is defined in Eq. (!26|) . 

To perform the MFA calculations one starts by solving Eqs. (1231) and (1321) with pi and cr 
given by Eqs. (HSj) and (1461) . As outlined above, in contrast to the BEC of atomic gases, in 
the triplon problem the chemical potential ix is the input parameter, whereas the densities 
are the output ones. Bearing this in mind, we rewrite the main Eqs. f l23|) and f l32|) as 



A = /i + 2t/ (a - pi) , (48) 

po = A/f/-cr. (49) 

The system of coupled equations, fH5|) . (H6|) and (l48l) . (l49l) has to be solved for given T and 
p to evaluate the triplon density 

p(T<Te) = po + Pi = ^^, (50) 

which is proportional to the measured sample magnetization density M . Note that by 
formally setting in all above formulas cr = 0, one arrives at the HFP approximation and 
particularly 

A = p-2f/pi, Po = A/f/. (51) 



■^Note that the second integral in Eq. (pij) can be treated with the Vehman formula, see e.g., H. Kleinert 
Hagen and V. Schulte-Frohlinde, Critical properties of ^''-theories, World Scientific Publishing Company 
(200f). 
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2.3. The critical temperature and triplon density 

Before discussing the normal T > phase we evaluate the temperature of BEC transi- 
tion. It is well-known that in the MFA the system of interacting Bose condensate particles 
at T —;■ Tc behaves like an ideal gas. In fact, assuming po{T — Tc) = 0, pi(T — )■ Tc) = p, 
a{T Tc) = 0, A(T T^) = 0, and Ek = Sk one concludes from Eqs.(ll5]),(ll8]) that 



/■ d^k 1 /i 

= ] (2vr)3exp(/3,e,)-l = 2Z7- ^^^^ 

With given p = /i^S'-f^ext — ^st and coupling constant f/, the critical temperature Tc = l//3c 
can be found as a solution of Eq. (152|) . Note that for the parabolic dispersion Sk = fcV2m, 



the integral (152|) can be evaluated analytically giving following well known relation [41 

2/3 



m \2C{3/2)U J ' 



(53) 



where ({x) is the Riemann function. 

As it has been underlined in Introduction, the bare dispersion of magnons, Sk, plays the 
crucial role in determining Tc. Figure H] presents the transition temperature as a function of 
magnetic field ifext and corresponding triplon density p, that is the sample magnetization, 
for parabolic Sk = k"^ /2m and "relativistic" 



Sk = JAI + J^k^A - Ast (54) 



bare dispersion, typical for a system with a gapped spectrum (23|, |42]. Here the effective 



exchange parameter J = 2^ iS.^^jm is chosen to match the parabolic and the relativistic 
Sk at small k. The dashed line obtained directly from Eq. fl53p displays T^ ~ p^/^ behavior. 
The solid one is a result of numerical solution of Eq.f l52p with the dispersion in Eq. fl54p and 
shows a crossover from Tc ~ j?!"^ at lower to Tc ~ (P'^ at higher temperatures in agreement 
with the experimental datajj The decrease in Agt enhances the role of relativistic features 
in the spectrum and leads to a faster deviation from the Tc ~ p^^^ behavior, as can be 
seen in Fig|2J Here and below we mainly use the set of input parameters as m = 0.0204 



K~ , Ast = 7.1 K, f/ = 313 K, and g = 2.06 [3lJ valid for the weakly anisotropic quantum 
antiferromagnet TlCuCls. 

In the normal phase the symmetry in Eq.(|3]) is not broken and the Bogoliubov shift is 
not needed either. Here the anomalous density vanishes, a{T > Tc) = 0, and hence, both 
approximations, HFB and HFP coincide. As a result we obtain for the triplon density 

d^k 1 

(2^)3 exp (/3wfc) - 1 (55) 



p{T >Tc) = I 



■*We mention that for linear dispersion Sk ~ k, the exponent would be equal to 1/3, and, therefore, in 
the experimental regime, the triplon gas is between the nonrelativistic and strongly relativistic realizations. 
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Figure 2: (Color online) The critical temperature of BEC of triplons Tc as a function of external field (left 
panel) and corresponding triplon density p = /i/2[/ (right panel) for two values of the singlet-triplet gap. 
To understand the role of the gap in the magnon spectrum, we present the results for Agt = 4 K, with the 
reduced exchange parameter yielding the same effective mass as for Agt = 7.1 K. The gap values are marked 
near the lines. The solid and dashed lines are for the relativistic and parabolic dispersions for Agt = 7.1 K, 
respectively. 



with Uk = Ek - fi + 2Up = Ek - yUeff- Similarly, by setting A = po = a = 0, pi = p, Ek = ujk 
in Eq. (l47|) one obtains following equation for the energy per unit cell 

The density of triplons as a function of temperature, evaluated in the HFP and HFB 
approximations is presented in Figj3l It is seen that the former leads to a discontinuity in 
the magnetization near the critical temperature, while the latter gives a continuous behavior 



in accordance with the experimental data [3jJ, |43|. However, in the condensate phase, the 



triplon density is higher in the HFB than in the HFP approximation. Therefore, the validity 
of the description of the BEC in the HFP approximation can be checked in the magnetization 
measurement experiments. 

3. Macroscopic properties: specific heat and magnetic susceptibility 

3.1. General expressions with nonzero anomalous averages 

Triplon contribution to the constant volume specific heat C^, can be calculated by dif- 
ferentiation of the energy with respect to the temperature at given chemical potential, that 
is at given external field 4J, |45| : 



a(i/ext,T) = ||, (57) 

where the energies for condensate and normal states are given by Eqs. (jUj) and (jSSD, 
respectively. 
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Figure 3: (Color online). Comparison of the HFB (solid lines) and the HFP (dashed lines) results for the 
triplon density. The HFB approach shows a continuous behavior, which fully agrees with the experimental 
data [31I, while the HFP approach leads to the discontinuity. The corresponding magnetic fields i/oxt 
are marked near the plots. 



Before discussing the magnetic susceptibility, we note that the macroscopic properties 
of the systems are related to their response to external fields. An example is given by the 
isothermal compressibility 



If — > 00 the system becomes unstable j45|, |46| since an infinitesimal fiuctuation of pressure 
P will lead to its collapse or explosion. As to the triplons their density is proportional to the 
magnetization M, while the i^ext simulates the pressure, and hence the relevant parameter, 
which determines the stability of the system with respect to the magnetic field, is the 
susceptibility 

dM 



The susceptibility can be calculated with the triplon density as: 

x{H,^„T) = {gfiBf^. (60) 

We will omit ifgxt from the arguments of x{Hext,T) and Cv{Hcxt,T) below. 

We begin with the normal phase, where po = = ^ = 0. The derivative of the density 
with respect to the temperature can be evaluated here with Eq.f lSSp as 

dp /3S, 

Calculating the derivative of Eq.( 15B]) with respect to T one obtains: 

a(r > T,) = -^3 + 2f/5i^- (62) 
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Here we introduced dimensionless quantities 

Si = Y,^kfBi^k), S2 = uY,fBi^k), Ss = (3Y,^lfBi^k), (63) 

k k k 

and used the relation (9/b(w)/5T = —(3ujf'^{uj), where (w) = dfB{uj)/duj = —(3exp{/3uj)f^ {u) 
and Uk is given in Eq. fl22|) . 

Similarly, taking the derivative of the self-consistency Eq. (l55|) . which we present here in 
the form: ^ 

^^^^ ^ ^ exp[/3 {e,-^ + 2Up{f,))]-l ^^^^ 

with respect to /i and solving the equation for dp/dfi one finds dp/dp = S2/U{2S2 — 1). As 
a result, we obtain with Eq. (l60|) the normal phase susceptibility: 

,(r>r,^(«^^. ,65) 

Now we proceed with the condensate phase, where the dependence of A and the cor- 
responding normal and anomalous densities on temperature and magnetic field should be 
taken into account. 

Here C„(T) is obtained by taking the derivative of (T < T^) given by Eq. (147|) . The 
latter can be rewritten as: 

E{T<T,) = E2-^ + U{pI~ 2p,a) , (66) 

where we used the relation po = fi/U — 2pi + a and present E2 defined in Eq. fl37Il as: 

IStt 



i^. = ^^^^^ + E^^-/-(^^)- (67) 



As a result: 

C.iT < T.) = ^ - 2Up,^ -2u(p- ^) ^. (68) 

We begin with calculation of dA/dT, which is the key ingredient in the specific heat. 
It is obtained by differentiating both sides of the main equation f HS]) with respect to the 
temperature, and for known dA/dT the other two derivatives dpi/dT and da/dT can be 
evaluated directly from Eqs.f l45p . (146|) . As a result one obtains the derivatives necessary to 
evaluate the specific heat: 

OA 2U(5 



dT A/Arr73/2 _ 



J2i^k + 2A)f'^{E,), (69) 



k 



dpi m\J Am dA 

W ^ 27r2 dT ' ^ [dT yEl" ^^"^ ' E| 



da 3my/ Am dA ^ ^ 
df ^ 27^ df ~ ^ 



dAfe.A^,^, , 4 

^ ( ^-F(E,) + 4rfB{E,) ] - /3AUE,] 



'k 



9^ f^^^T(F)^'' 



(70) 
(71) 
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where we introduced notation for the frequently used expression: 



T{E,)^l^ + aE,). (72) 



The expression for 

dE2 A{Amf^^dA sr^fdA 



dT 37r2 dT 



J2 (1^^^-^ (^^0 - ^ElaE,)^ (73) 



completes the set of derivatives necessary for evaluation of the specific heat. 

The derivative dA/dfi can be found by differentiating both sides of the main equation 
( HSj) with respect to /i. As a result, 

^ = I (74) 

As to dp/dfj,, needed for evaluation of the magnetic susceptibility in Eq.( l60|) . it is obtained 
by differentiating the equation p = (A + /i) /2U, with respect to /i. This yields 

Below we apply these equations to the cusps in the specific heat and susceptibility and to 
the qualitative effects such as the instability in strong magnetic fields. We will show that 
at a given temperature, there exists a critical field H^l^^{T), such that when i/cxt approaches 
> ifg^j(T), the magnetic susceptibility, x{Hcxt,T) diverges. 

3.2. Cusp in the specific heat and magnetic susceptibility nearT^ 
The cusp in the specific heat defined as 

AC,, = lim CJT) - lim CJT) (76) 

is an interesting quantity in the theory of phase transitions. In accordance with the Ehrenfest 
classification, a phase transition with the discontinuity in C„ near the transition point, is 



the second order one. Particularly, it is well-known that |44l . |45| for the ideal gas AC^, near 



the transition into BEC is zero i.e. the specific heat is continuous. We shall illustrate this 

pecific heat per particle of the i 

^^-^^'-^CWWyTj ' (77) 



fact for completeness. The specific heat per particle of the ideal Bose gas is given by |44j . |45 

3/2 



c^°(r > r.) = 



4|3/2(^) 



9g3/2(^ 
4^1/2 



(78) 



where gp{z) is defined as: 



^P-i 



-) = r(rtX '"F^^- '''' 



related to the Riemann function as ({p) = gp{z = 1). When T — Tc + 0, the fugacity 
z = exp(Bii) < 1 tends to unity, i. e. lim z = 1 and the second term in flTHl) vanishes, 

T-i-Tc+O 

because of the divergence in gi/2, i.e. gi/2{z) ~ (1 — 2;)"^/^ in this region, while the first 
term exactly coincides with fl77p . As a result, 

AC^^=o = 0, (80) 

and hence the specific heat of the ideal gas is continuous although being plotted as a function 
of temperature C^^°(T) behaves similarly to the A curve. 

We proceed with calculation of AC^,. Bearing in mind that for T > Tc HFB and HFP 
approximations coincide, from f l62|) we obtain 



lim C^,{T) = -PcY.elUe,). (81) 

k 

For T < Tc, is given by Eqs. f l68|) -f l72|) . Assuming in the last equations Ek = Sk and 
A = 0, one finds 

9A 2PcU 
hm — = — > efc efc, 82 

T-,T.~0dT l + 2[/^^(£fc) V 

k 

which is finite. Using ([82]) in ([70])-([72]) and setting in ([68]) p = /i/2?7, po = 0, pi = p, cr = 0, 
one obtains for the maximum value of the specific heat max{C^™(T)} = \\mT^Tc-oCv{T): 



M) 



^hm__^a(T) = 5^(/i + 5,) U{e,)^-p,e,f^{e,)\ . 

k ^ 

Subtracting (IHTl) from (l83l) we finally obtain 

AC„ = 5^ f (^fc) (/i + ^fc) 1^ - /ute/fil^fc) j , 

k ^ 

where dA/dT is given by Eq. (182|) . 

In a quite similar way one can calculate the cusp in the susceptibility: 

. {gf^B^ 1 
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3.3. Parabolic dispersion, Sk = k'^/2m 

For the parabolic dispersion, calculations can be done analytically, as presented below. 
The specific heat is expressed as: 

dA 

max{a,(T)} = -(3,Ri - ^P,R2 + ^ {R2 + fiR-s + R5) , (86) 
lim a(T) = (87) 

dA 

AC, = max{a(r)} + /3,Ri = -^{R2 + /ii?3 + R^) - /i/3ci?2, (88) 
and the cusp in the susceptibility has the form: 

The derivatives dA / dT and dA / dfi can be simplified to: 

OA _ 2PJJR2 
df ~ 1 + 2UR3 ' 
dA 1 



(9/i 1 + 2UR3 
In these formulas we used following notations: 



(90) 
(91) 



Ri = -/3.^/|(..)4e^--- = - ^^^/"^^jJ/'^ T.^, (92) 

k 

D Q sr^ f2 \^ B^sk C(3/2)r (5/2) 3 

R2 = -f)c2_^fB{^k)£ke'^'' = Kt, (93) 

k 

= E(/^(^^) + /-(^^)/^^-) = -^^^^f, (94) 

k ic 

= 5^/B(e.)e. = ^^^^^^|^T.4, (95) 



where we introduced kx = VmTc for the characteristic thermal wavevector of triplon at the 
transition temperature. Taking into account that R5 is the triplon concentration, we find 
that Eq. fl96|) is equivalent to Eq. fl53|) . 



4. High- field instability 



It is well-known that |46|, |47[ the dynamic stability of equilibrium system is determined by 
its excitation spectrum: When the latter becomes imaginary the time evolution of excitations 
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Figure 4: The LHS of Eq. (|10ip for different values of 77. The real solutions of this equation correspond to 
the intersection of F{Z, -q) with the Z = axis. 



changes qualitatively. In the condensate state the spectrum is given by = ^/^\/ + 2 A 
where A is the solution of nonlinear algebraic Eq. fl48p which may be rewritten as follows 

4C/(Am)3/2 I 



A = n + 



37r2 



vr 



gfc + 2A 
Eh 



(97) 



Obviously, for some set of system parameters U, m and external fields ifcxt, Eq.( l97|) may 
have no real solution and hence the spectrum, as well as the sound speed, become imaginary. 
For example, in our previous work we have shown that for parameters 3l| this may happen, 
at H^^^ = 12.5 T at T = with the increasing in H^^^ with T. Below we investigate 
the magnetic susceptibility near the instability point and show that it goes to infinity, as 
expected in these cases, providing the direct experimental test for the instability. This effect 

/icr the denominator of Eq. 
the main equation f p7|) has 



can be seen from the fact that when ifext = H^it^ i- ^- 



f in|) becomes zero, and hence dA/dfi diverges. For > 
no positive solution and hence the sound speed c = -\/A/ 
illustrate this analytically for T = 0. 

For zero temperature the derivative dA/dfi in Eq. 

dA _ 1 
diJ. ~ 2C/\/Am3/2 

where A is the solution to the following equation: 

4C/(Am)3/2 



m becomes complex. Below we 



31) is written as 



(98) 



A = /i + 



Introducing dimensionless variables rj 
equations as: 

dA _ 



37r2 
and Z 



(99) 



A//i one can rewrite the last two 



(100) 



4Z3/2y^ 
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37r2 



-1 = 0, 



(101) 



where t] is an input parameter. In FigJH 77) is plotted as a function of Z for various 
values of 77. It is seen that for 77 larger than some critical the second term in fllOip 
dominates, F{Z, 77) is always negative, and hence there is no real solutions. If the maximum 
value of F{Z, rj), (maxF(Z, rj)), for a given 77 is negative the Eq. fllOip . has no real solution. 
Otherwise, when max(F(Z, 77) is positive, F{Z, rf) intersects the Z-axis and the equation has 
at least one positive solution. Thus, the boundaries and Z^r are defined by the coupled 
equations: 

^^3/2 — 

r/er) = Zcr ^-P^ -1 = 0, (102) 



OF 
dZ 



= l-hl?^ = 0, (103) 

TV 

Z=ZcT,r]=ilci 



giving 7/cr = vr^/12 = 8.1174 and Zcv = 3. By comparing (11031) and (llOOl) one concludes that 
when 7] approaches rj^r, dA/d^ goes to infinity and so does the magnetic susceptibility given 
by Eq. (175]) i.e. x can diverge. When rj exceeds 7/cr, solutions of Eqs. (lUUi) . (llOip acquire an 
imaginary part. Bearing in mind that rj = fim^U'^ = {fiBgH^^t — Ast)'m^U'^, one concludes 
that even at T = 0, if the Hcxt is strong enough the speed of sound c = A/m = a/ j^Z/m 
becomes complex and the Bose condensed system of triplons displays dynamical instability. 
By expanding Eq. OlOip in the vicinity of the point Z = Zcr,7/ = rjcr, we obtain in this 
region for the real and imaginary parts of the speed of sound: Re(c) = ^J^jj^^Jm and 
Im(c)/Re(c) = — — ficrUm^^'^ /tt"^, where = Vcr/i'm^U'^)- 

Similarly it can be shown that for any given T < there is a maximal value of the 
external magnetic field such that the system posses dynamical instability at this point 
i.e. x(-f^ext)^") ~^ ^ ^'^d simultaneously the sound speed acquires an imaginary part. 
This fact is illustrated in FigsOandO In Fig|S]the susceptibility x{Hext,T) is plotted as a 
function of -ffext- In FigElthe solution of the main equation, more precisely, the components 
of the sound speed c are plotted as a function of iJcxt for several temperatures. It is clearly 
seen that for a given temperature the x{Hext,T) diverges and the sound speed becomes 
complex at the same Hext- 

Note that in the HFP approximation Eq. OlOip being written as 

FnFp{Z, r]) = Z+ -1 = (104) 

has a positive solution for any positive 7/. Thus this approximation precludes the appearance 
of the instability. 

The experiment related results from the observed instability are presented in Figj5] and 
Figjni Figure |5] shows the divergence of the susceptibility. Figure |6] shows how the imaginary 
part in the velocity appears and the real part of the velocity changes if the field becomes 
stronger than the threshold value. This is a qualitative effect, which should demonstrate 
itself in the experiment. Here an important comment on the values of the fields where 
the singular behavior of A can be observed is in order. A comparison of the calculations 



presented in Fig|3l experimental results [311 and theory taking into account the exact 
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Figure 5: The magnetic susceptibility xC^y -ffext) for the temperatures marked near the plots. When external 
magnetic field reaches the critical value -ffext susceptibility diverges. 



spectrum of triplons 28|], shows that our model spectrum leads to an overestimate of the 
tiplon concentration by about a factor of 1.5. Therefore, it underestimates the threshold 
field, leading to the result that fields of the order of 20 T will be necessary 48|] to put the 
condensate in the unstable part of the phase diagram. 

We mention two characteristic features of the instability considered in this Section. First, 
it is not a hydrodynamical instability arising in the inhomogeneous flow of the Bose-Einstein 



condensate [49|-|52|. Second, it is not directly related to the phonon-phonon scattering [53 



since we are still in the mean field approximation, and the Hamiltonian we consider does 
not contain higher-order products of the phonon operators. Nevertheless, both these effects 
can be studied taking into account the anomalous averages in the theory. 

Summarizing this Section we conclude that there are two characteristic values of H^xt- 
First, when the triplons form the BEC, Hf^'^ and second, if^^^ > H^J^'^, when the BEC 
displays dynamical instability. The phase diagrams for different repulsion parameters U for 
relativistic triplon dispersion are shown in FigJTl The localization of the stable BEC phases 
on the (T, i^cxt) plane depends on the model parameters (m, U, Ast). One can conclude that 
increasing of U decreases Hf^^ and makes the BEC less stable, as expected from the fact 
that the instability is caused by the magnon-magnon repulsion. 



5. Conclusions 

We investigated macroscopic properties of the triplon condensates taking into account the 
anomalous density averages. We show that these averages have an important quantitative 
effect on the specific heat and qualitatively modify the magnetic susceptibility of these 
systems. The qualitative result of our consideration is related to the dynamical instabilities 
arising in the condensate. These instabilities are seen as the divergence in the magnetic 
susceptibility and changes in the sound of the Bogoliubov mode when the external magnetic 
field becomes stronger than the corresponding temperature-dependent critical value. In this 
domain, the speed of sound acquires an imaginary part and dependence of its real part on 
the external field becomes weak, almost fiat. These two predictions can easily be verified 
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10 12 14 16 18 20 



Figure 6: The components of the sound speed c for different temperatures, real for Hcxt < ^^oxti ^^'^ complex 
out of this range. Real parts Re(c) (solid lines) are positive and Im(c) (dashed lines) are negative. Near the 
threshold Im(c) demonstrates to the mean- field square root behavior as discussed below Ea. (|103p . 




Figure 7: The instability borders (dashed lines) as well as the condensate formation borders (solid lines) for 
different U. With the increase in U the area corresponding to the stable condensate decreases. 



experimentally in laboratory studies of the static magnetic susceptibility and in the neutron 
scattering measurements in various magnetic fields. 



To conclude this paper, one general comment should be given. Several authors [5J, |55 



casted doubts on the applicability of the Bose-Einstein condensation approach to magnetic 
transitions. The irrefutable proof of the validity of this approach can be given by the 
observation of the Josephson effect, as it was done for magnons in the superfluid He^ ^ 
and although recently proposed for magnetic insulators 56[, not yet observed for the triplon 
realization. However, we believe that observation of the effects predicted in this paper which 
are solely based on the improved theory of the Bose-Einstein condensate, will resolve this 
controversy. 
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